ORBITAL LINEARIZATION OF SMOOTH COMPLETELY 
INTEGRABLE VECTOR FIELDS 



NGUYEN TIEN ZUNG 

Abstract. The main purpose of this paper is to prove the smooth 
local orbital linearization theorem for smooth vector fields which admit 
a complete set of first integrals near a nondegenerate singular point. 
The main tools used in the proof of this theorem are the formal orbital 
linearization theorem for formal integrable vector fields, the blowing-up 
method, and the Sternberg-Chen isomorphism theorem for formally- 
equivalent smooth hyperbolic vector fields. 



1. Introduction 

The main purpose of this paper is to show the fohowing orbital lineariza- 
tion theorem for smooth (C°°) vector fields which admit a complete set of 
first integrals near a nondegenerate singular point: 

Theorem 1.1. Let X be a smooth vector field in a neighborhood of O = 
(0, ... ,0) in M", which vanishes at O and satisfies the following conditions: 
i) (Complete integrability) : X admits n — 1 functionally independent smooth 
first integrals Fi, . . . , i.e. X{Fi) = . . . = = and dFi A ... A 

dFn-i ^ almost everywhere. 

a) (Nondegeneracy 1): The semisimple part of the linear part of X at O is 
non-zero, and the oo-jets of Fi, . . . ,-Fn-i o-t O are funtionally independent. 
Hi) (Nondegeneracy 2): If moreover is an eigenvalue of X at O with 
multiplicity k > 1, then the differentials of the functions Fi,...,Fk are 
linearly independent at O: dFi{0) A ... A dFk{0) / 0. 
Then there exists a local smooth coordinate system (xi, . . . in which X 
can be written as 

(1.1) X = FX^'^\ 

where X^^^ is a semisimple linear vector field in (xi,...,x„), and F is a 
smooth first integral of X^^\ i.e. X^^\F) = 0, with F{0) = 1. 

The above theorem is in fact more than mere orbital linearization: not 
only that X is orbitally equivalent to its linear part X'^^\ but also the factor 
F in the expression in a normalized coordinate system is a first 
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integral of X and X'^^\ In this kind of linearization is called geometrical 
linearization. 

The formal and analytic case of the above theorem also holds and was 
shown in [9] in a more general context of integrable non-Hamiltonian sys- 
tems of type {p,q), i.e. with p commuting vector fields and q common first 
integrals, where p + q = n is the dimension of the manifold. The vector 
fields that we study in this paper are integrable of type (l,n — 1), i.e. just 
one vector field and n — 1 first integrals. 

The nondegeneracy condition in Theorem 11.11 is a bit stronger than the 
nondegeneracy condition in [9] : in [9] the (formal or analytic) vector field X 
is called integrable nondegenerate if it satisifes the above conditions i) and 
ii), without the need of condition iii). In fact, in the formal and analytic 
case, condition iii) is a simple consequence of the first two conditions and 
the theorem about the existence of (formal or analytic) Poincare-Dulac nor- 
malization [8l[9]. However, in the smooth case, we don't have a proof of the 
fact that condition iii) follows from conditions i) and ii) in general, though 
we do have a proof of this fact for dimension 2. 

The rest of this paper is organized as follows. Section 2 is devoted to 
some preliminary results, including the classification of the nondegenerate 
singularities of completely integrable vector fields into (strong/weak) ellip- 
tic and hyperbolic cases (Lemma 12. 2p . and the normalization up to a fiat 
term (Proposition 12. 3p . These preliminary results are used in the proof of 
Theorem 11.11 which is presented in Section 3. Finally, in Section 4, we show 
that, at least in the case n = 2, condition iii) Theorem 11.11 is a consequence 
of the first two conditions, and can be dropped from the formulation of the 
theorem (Theorem 14. ID . We conjecture that condition iii) is redundant in 
the higher-dimensional case as well. 

This paper is part of our program of systematic study of the geometry 
and topology of integrable non-Hamiltonian systems. In particular, Theorem 
14. H which is a refinement of Theorem 11.11 in the case of dimension 2, is the 
starting point of our joint work with Nguyen Van Minh on the local and 
global smooth invariants of integrable dynamical systems on 2-dimensional 
surfaces [lOj . 



2. Preliminary results 

2.1. Adapted first integrals. We have the following simple lemma, which 
is similar to the well-known Ziglin's lemma [7j 

Lemma 2.1. Let Gi, . . . ,Gm be m formal series in n variables which are 
functionally independent. Then there exists m polynomial functions of m 
variables Pi, . . . , Pm such that the homogeneous (i.e. lowest degree) parts 
of the formal series of Pi(Gi, . . . , Gm), ■ ■ ■ , Pi{Gi, . . . , Gm) are functionally 
independent. 
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The proof of the above lemma follows exactly the same lines as the proof 
of Ziglin of his lemma in |7] , and our situation is simpler than the situation 
of meromorphic functions considered by Ziglin. 

Let X be a smooth completely integrable vector field with a singularity at 
O. We will say that the smooth first integrals Fi, . . . , Fn-i of X are adapted 
first integrals if 

(2.1) dHiA...AdHn-ij^O a.e., 

where Hi = F^^'''^ denotes the homogeneous part (consisting of non-constant 
terms of lowest degree in the Taylor expansion) of Fi at O. Using the 
above lemma to replace the first integrals Fi, . . . , Fn-i of X by appropriate 
polynomial functions of them if necessary, from now on we can assume that 
Fi, . . . , Fn-i are adapted. 

2.2. The eigenvalues of X. The fact that X admits n — 1 first integrals 
implies that the X is very resonant at O. More precisely, we have: 

Lemma 2.2. Let (X, Fi, . . . , F„_i) he smooth nondegenerate at O, i.e. they 
satisfy the conditions of Theorem li.il Then the linear part of X at O is 
semisimple, and there is a positive number A > such that either all the 
eigenvalues of X at O belong to XL, or all of them belong to y/—lXL. 

Proof. We can assume that Hi, . . . , Hn~i are functionally independent, where 
Hi denotes the homogeneous part of Fj. The equality X{Fi) = implies 
that 

(2.2) X''{H,) = X« (i/i) = Vi = 1, . . . , n - 1, 

where X^^^ is the linear part of X, and is the semisimple part of X^^) 
in the Jordan-Dunford decomposition. We can write 

(2.3) ^^^ = EV.9^ 

1=1 * 

in a complex coordinate system. Recall that the ring of polynomial first 
integrals of X^ILi ^i'^iW' ^® generated by the monomial functions XXi= 
which satisfies the resonance relation 



(2.4) ^a,Ai = 0. 

2 = 1 

The fact that H\,...,Hn-\ are independent implies that Equation (j2.4p 
has n — 1 linearly independent solutions which belong to , which in turn 
implies that there is a complex number A such that Ai, . . . , A„ G XL. Remark 
that if the spectrum of X^^ contains a complex eigenvalue Ai G C \ (M U 
\/— IM), then its complex conjugate Ai is also in the spectrum because X is 
real, and Ai and Ai cannot belong to AZ at the same time for any A. Thus 
any eigenvalue of X^'^ is either real or pure imaginary. If there is one real 
non-zero eigenvalue, then we can choose A G R_|_, otherwise we can choose 
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A G \/— Notice that A / because at least one eigenvalue of X^^ is 
non-zero by our assumptions. 

The common level sets of Hi, . . . are 1-dimensional almost every- 

where, and since both X^-^'^ and X^'^ are tangent to these common level sets, 
we have that X^^-* A X'^^ = 0, which implies that X^^^ is semisimple, i.e. 

With the above lemma, we can divide the problem into 4 cases (here 
R* =M\{0}): 

I. Strongly hyperbolic (or hyperbolic without eigenvalue 0): Aj G AM* Vi. 

II. Weakly hyerbolic (or hyperbolic with eigenvalue 0): Aj G AM* Vi > 
A; > 1, Ai = . . . = Afc = 0. 

III. Strongly elliptic (or elliptic without eigenvalue 0): Aj G -v/^^AM* Vi. 

IV. Weakly elliptic (or elliptic with eigenvalue 0): Aj G \/— 1AM* Mi > k> 
1, Ai = . . . = Afc = 0. 

2.3. Linearization up to a flat term. Using the geometric linearization 
theorem of |9j in the formal case, we get the following proposition: 

Proposition 2.3 (Linearization up to a flat term). Assume that X satisfies 
the hypotheses of Theorem Then there is a local smooth coordinate 

system (xi, . . . ,2;„) in which X can he written as 

(2.5) X = FX(^) + flat 

where X^^^ is the linear part of X in the coordinate system (xi, . . . ,x„), F 
is a smooth first integral of X^^\ and flat means a smooth term which is 
flat at O. 

Proof. Denote by X (resp. Fi) the oo-jet of X (resp. Fi) at O: X is a 
formal vector field (resp. function) at O. If (X, Fi, . . . , F„_i) is smooth 
nondegenerate at O, then {X,Fi, . . . ,F„_i) is a nondegenerate formal inte- 
grable system of type (1, n— 1) at p. According to the geometric linearization 
theorem of [9] , this formal integrable system can be linearized geometrically, 
i.e. there is a formal coordinate system (xi, . . . , x^) in which we have 

(2.6) X = FX^^^ 

where X^^^ is the linear part of X in the formal coordinate system (xi, . . . , x„), 
and F is a formal first integral of X^^\ By the classical Hilbert-Weyl theo- 
rem (see, e.g., [3] ), we can write 

(2.7) F = /(Qi(xi, . . . , x„), . . . , Qm{xi, ■■■ , x„)) 

where / is a formal series and Qi{xi, . . . , x„), . . . , Qmixi, • • • , x„) are ho- 
mogeneous polynoms generate the ring of polynomial first integrals of X^^\ 
Using Borel theorem, we get a smooth coordinate system (xi, . . . , x^) whose 
c«-jet is (xi, . . . ,Xn), and a smooth function / of m variables whose oo-jet 
is /. Put 

(2.8) F(xi, . . . , Xn) = fiQlixi, ... , Xn), ■■■ , Qm{xi, Xn))- 
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Then Equations ([M]), dZZ]) and imply that X = FX^^ + flat in the 
smooth coordinate system {xi, . . . , Xn)- □ 

2.4. Reduction to the case without eigenvalue 0. Assume that X has 
zero eigenvalue at O with multiplicity k, and dFiA. . .AF^ = 0, i.e. we can use 
Fi,. . . ,Fk as the first k coordinates in our local coordinate systems. Since 
the vector field X preserves xi, . . . we can view it as a /c-dmiensional 
family of vector fields on (n — A;)-dimensional spaces 

(2.9) Uc,,...,c, = {Fi=ci,...,Fk = Ck} 

(for ci , . . . , Cfc small enough) . It follows from the usual implicit function 
theorem that on each Uci,...,Ck there is a unique point Oci,...,Cfc such that 
^{Oci,...,Ck) = 0) and moreover the point Oci,...,Cfe depends smoothly on 
ci, . . . , Cfc, the eigenvalues of X at ci, . . . ,Ck, are non-zero. It also follows 
from the formal independence oi Fi, . . . , Fn at O, that the functions func- 
tions -Ffc+i, . . . ,Fn are formally independent at every point Oci,...,Cfc provided 
that ci , . . . , Cfc are sufficiently small. In other words, we have a /c-dimensional 
family of nondegenerate singularities of smooth completely integrable (n—k)- 
dimensional vector fields -'^^ci,...,Cfc- In order to normalize X, it suffices to 
normalize Xcj,...,Cfe in a way which depends smoothly on the parameter. 

3. Proof of Theorem 11.11 

We will always assume that the vector field X satisfies the hypotheses of 
Theorem ll.il The fact that the linear part of X is semisimple is established 
by Lemma 12.21 In view of Subsection 12. 4^ it suffices to prove Theorem 11.11 
for the cases without zero eigenvalue, by a proof whose parametrized version 
also works the same. 

3.1. The hyperbolic case. Assume that X is hyperbolic without eigen- 
value 0. According to Proposition 12.31 we can write X = Y + flat, where 
Y = FX^^^ is a smooth hyperbolic integrable vector field in normal form. 
Since X and Y are hyperbolic and coincide up to a flat term, Sternberg- 
Chern theorem [Bl [2] says that X is locally smoothly isomorphic to Y, 
i.e. there is a smooth coordinate system in which X can be written as 
X = where F is a smooth first integral of X^-^\ Theorem 11.11 is 
proved in the hyperbolic case without eigenvalue 0. 

3.2. The elliptic case. In this subsection, we will assume that all the 
eigenvalues of X at O are non-zero pure imaginary. Using Proposition 12. 3( 
we can assume that X = + flat in a local smooth coordinate system 
(xi, . . . ,Xn), where F is a smooth function such that F{0) = 1. Put Y = 
X/F. Then Y has the same first integrals as X, and 

(3.1) Y = X(^) + flat. 

The fact that X is of strong elliptic type implies immediately that the 
dimension n is even, the eigenvalues of X at O are it-v/— loi, • • • , ±\/— lOn/2 
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where oi, . . . , a„/2 are positive real numbers, and we can choose the coordi- 
nates (xi, . . . , Xn) such that 

(3.2) X« = Va,(x2._i^^ X2^J. ). 

^ dx2i dx2i-l 



According to Lemma 12.21 we can choose A > such that ai/A, . . . ,an/2/A 
are natural numbers whose greatest common divisor is 1. 

Lemma 3.1. Locally near O all the orbits of Y = X/F (except the fixed 
point O) are periodic, with periods which are unifromly bounded above and 
below. 

Proof. The vector field {dFi A ... A dFn-i)j{^ A ... A ^) is tangent to 
Y, and therefore it is divisible by Y, i.e. we can write 

d d 

(3.3) (dFi A . . . A (ii^„_i)j(— A ... A ^) = GY, 

where G is a smoth non-flat function at O. Notice that the singular locus of 
the map (Fi, . . . ,Fn-i) ■ U — )• M""^, where U 3 O is a small neighborhood 
of O in R", coincides with the zero locus of G. 

It is clear that, by continuity, the set of all points x £ U such that the 
orbit of Y through x is periodic of period < Svr/A is a closed subset of U. 
We want to show that this set is actually equal to U (provided that U is 
small enough). Consider the singular locus 

(3.4) S = {xeU \ G{x) =0} = {x eU \ dFi A... A dFn-i{x) = 0} 

Since G is non-flat at O, we can choose a coordinate system (zi, . . . ,z„) 
which is a linear transformation of the coordinate system (xi, . . . , x„), such 
that the homogeneous part G^'^^ of G has the form 

(3.5) G^^) = z^ + ... 

which means that 7^ in ?7. Because does not vanish in U , by 

the classical Rolle's theorem on each line {z2 = const, . . . ,Zn = const} in 
U there are at most k zeros of the function G, the singular locus S is of 
dimension at most n — 1, and the function G is not flat at any point of S. 

Take a point = {zi = e, Z2 = 0, . . . , z„ =0) G U with e > small 
enough. Take the (n — l)-dimensional ball B"'~^{q^, e^) of radius e^' which 
is orthogonal to the vector Y{q^) at the point q^ in the coordinate system 
(zi, . . . , Zn), for a certain positive number K to be chosen below. Denote 
by 4> the Poincare map of the flow of y on D. A-priori this map does not 
necessarily fixes the point g^. But due to the fact that Y = X^^^ -|- flat and 
the flow of X^^^ is periodic, the distance from q^ to </>((?e) is smaller than 
e^"*"^ (provided that e is small enough). We can choose K large enough so 
that the restriction of the map (Fi, . . . ,F„_i) to B"'~^{q^,e^) is injective. 
Due to the invariance of the functions with respect to the vector field 
Y, the Poincare map (p, we also have the points q^ and 4>{qe) have the same 
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image under the map (Fi, . . . , But this map is injective on the baU 

e^) which contains these two point, so in fact these two points must 
coincide, i.e. wc have = (f>{qe), and the orbit of the flow of Y through the 
point Qe is a periodic orbit, and the period of this orbit is equal to 27r/A plus 
a small error term which tends to faster than anu power of e when e tends 
to 0. 

Denote by V the path-connected component of f7 \ 5" which contains the 
points Then the orbit of Y through any point q & V is also periodic 
and its period is close to 27r/A (the difference between the period and 27r/A 
tends to uniformly when the radius of U tends to 0). This fact can be 
proved easily by showing that the set of points of V which satisfies the 
mentioned property is closed and open in V at the same time: closed due 
to the continuity, and open because {Fi, . . . is regular in V and is 

preserved by the flow of Y. 

Let q E S he a point in the locus S which also lies on the boundary of V. 
Then by continuity, there is also a number T near 27r/A such that the time-T 
flow of Y flxcs the point q. In other words, the orbit of Y through q is also 
periodic, and the period is equal to T or a fraction T/moiT for some natural 
number m. As before, consider a (n— l)-dimensional ball 6) which is 

centered at q and orthogonal to Y{q), for some S > small enough. Consider 
the Poincare map (p oiY on B^^^{q, 6) corresponding to the time T (i.e. if 
the period of the orbit through q is T/m then consider the m-time itaration 
of the usual Poincare map). Since the intersection of B"'''^{q,S) with V 
contains an open subset of B^~^{q,S) whose closure contains q, and the 
Poincare map is identity on that open subset by the above considerations, 
the Poincare map on B"'~^{q,6) is equal to the identity map plus a flat 
term at q. On the other hand, this Poincare map must preserve the map 
{Fi, . . . , Fn-i)\B"-i(q,5), and the determinant of the differential of this map 
is not flat at q. It implies that the Poincare map must be identity in a small 
neighborhood of q in B^~^{q,S). Thus, we can "engulf" the set of points 
shown to have periodic orbits from y to a larger open subset of U which 
contains the boundary of V. Continuing this engulfing process, we get that 
the set of points in U having periodic orbits is actually the whole U. □ 

We will linearize Y = X/F orbitally, and then deduce the normalization 
of X from this linearization. In order to do that, let us consider the blow-up 
of R" at O, which will be denoted by 

(3.6) piE^U, 

where U 3 O is a neigborhood of O in W and p~^{0) = MP""^ is the 
exceptional divisor of the blow-up in E. We will need the following simple 
lemma, whose proof is steaightforward: 

Lemma 3.2. With the above notations, a function G or a vector field Z is 
flat at O in U if and only if its pull-back to E via the projection map p is 
flat along p~^{0) in E. 
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Denote by G (resp. Z) the pull-back of a function G (resp. vector field 
Z) via the projection map n : E ^ U of the blow-up. Then we have 

(3.7) Y = X(^) + Z 

in E, where Z is vector field which is flat along p~^(0), and X^^^ is a smooth 
periodic vector field in E of period 2tt/X. By Lemma \3.1\ the orbits of Y 
are closed, with periods close to the period of X^^\ Due to the flatness of 
Z along p^^{0), the period of Y at the points in E is equal to 27r/A plus 
a smooth function on which is flat along p^^{0). Projecting Y back to 
U and using Lemma 13.21 we get a smooth period function P = 27r/A -|- flat 
(which is invarant on the orbits) such that PY is periodic of period 1. In 
other words, PY generates a smooth T^-action. Using the classical Cartan- 
Bochner smooth linearization theorem for compact group actions, we find 
a smooth coordinate system, which we will denote again by (xi, . . . , x„), in 
which PX/F = PY is a linear vector field, i.e. in which we have 

(3.8) X = GX(^) 

where G is a smooth function and X^^^ is a linear vector field which satisfies 
Formula 

A-priori, the function F given by Proposition 12.31 is not a first integral of 
X (though it is a first integral of the linear part of X in some coordinate 
system), and so the function G = 2ttF/P\ in Formula (j3.8p is not a first 
integral of X in either. But we can normalize further in order to change 
G into a first integral. Indeed, by the arguments presented above, we can 
assume that G is a smooth first integral of X plus a fiat term, or we can 
write G = Gi(l -|- flat), where Gi is a first integral of X. Normalizing 
the new vector field Y = X/Gi instead of the old Y = X/F, we get a new 
smooth coordinate system in which PY = {27r/ X)X^^^ , where P is the period 
function of the new vector field Y, and it is a smooth first integral of the 
type constant + flat. In this new coordinate system we have that X is equal 
to its linear part times a first integral, and Theorem 11.11 is proved in the 
elliptic case i.e. without eigenvalue 0. 

Since our proof for the strong hyperbolic case and the strong elliptic case 
also works for smooth families of integrable vector fields. Theorem 11.11 is 
proved. 

Remark 3.3. According to a theorem of Schwarz [5], the smooth first integral 
F in the normal form in the elliptic case can also be written as 

(3.9) F = fiQiixi, . . . ,x„), . . . ,(5m(xi, . . . ,x„)), 

where Qi(xi, . . . , a;„), . . . , Qmixi, ■ ■ ■ , Xn) are homogeneous polynoms which 
generatedthe ring of polynomial first integrals of the linear vector field X^^^ . 
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4. The case of dimension 2 

The aim of this section is to show that condition iii) in Theorem 11.11 is 
redundant at least in the case of dimension 2. More precisely, we have: 

Theorem 4.1. Let X he a smooth vector field in a neighborhood of O = 
(0,0) in M^, which vanishes at O and satisfies the following conditions: 
i) (Complete integrahility): X admits a smooth first integral Fi. 
a) (Nondegeneracy): The semisimple part of the linear part of X at O is 
non-zero, and the co-jet of Fi at O is non-constant. 

Then there exists a local smooth coordinate system (x, y) in which X can he 
written as 

(4.1) X = FX^^\ 

where X^^^ is a semisimple linear vector field in {x,y), and F is a smooth 
first integral of X^^\ 

Proof. Remark that, in the case of dimension 2, there are only 3 possibilities: 
elliptic without zero eigenvalue, hyperbolic without zero eigenvalue, and 
hyperbolic with zero eigenvalue. The first two possibilities are covered by 
Theorem 11.11 It remains to prove Theorem 14.11 for the case when X has one 
eigenvalue equal to 0. By Proposition 12.31 we can assume that 

o o o 

(4.2) X = Fiy)x— + flati— + flat2j- 

ox ox ay 

in a smooth coordinate system {x,y), where flati and flat2 are two flat 
functions, and F{0) ^ 0. 
Denote by 

(4.3) S = {qeU\ X{q) = 0} 

the singular locus of X near O, where U denotes a small neighborhood of O 
in 'S?. The main point is to prove that S" is a smooth curve. If .S is a smooth 
curve, then we can write S = {y = 0}, the vector field X is divisible by x, i.e. 
Y = X/x is still a smooth vector field, which is non-zero at O, and therefore 
locally rectifiable and admits a first integral G such that dG{0) ^ 0. But G 
is also a first integral of X, so condition iii) of Theorem II. II is also satisfied, 
and Theorem 14.11 is reduced to a particular case of Theorem 11.11 
Denote by 

(4.4) Si = {{x,y)GU\ F{y)x + flati{x, y) = 0} 

the set of points where the ^-component of X vanishes. It is clear that 
S C Si, and Si is a smooth curve tangent to the line {x = 0} at O by the 
inverse function theorem. We will show that S = Si. 
Consider the cone 

(4.5) C = {{x,y)eU\\x\<\y\}. 
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Clearly, Si d C (provided that U is small enough). The non-flat first integral 
Fi of X in the coordinate system (x, y) has the type 

(4.6) Fi = f{y) + flat 

where f{y) = a^y^ + h.o.t. is a non-flat smooth function. It implis that the 
level sets of Fi in the cone C are smooth curves which are nearly tangent 
to the lines {y = const}. In particular, each level set of Fi in C intersects 
with Si at exactly 1 point. Since X is tangent to these level sets, and the 
^-component of X vanishes at the intersection points of these level sets 
with 5i, it follows that X itself vanishes at these intersetion points. But 
every point of is an intersection point of Si with a level set of Fi. Thus 
X vanishes on Si, and we have S = Si. □ 

Remark 4.2. Two-dimensional elliptic-like vector fields, i.e. those vector 
fields whose orbits near a singular point are closed, are also called centers 
in the literature. There is a recent interesting theorem of Maksymenko [1] 
about the orbital linearization of the center, without the assumption on the 
existence of a first integral, but with an assumption on the periods of the 
periodic orbits. Maksymenko's theorem is similar to and a bit stronger than 
the elliptic case of Theorem [4T] because his assumptions are weaker, and the 
conclusions are the same. His proof is also based on the formal normalization 
and the blowing-up method. 

Remark 4.3. Some of the arguments of the proof of Theorem 14.11 are still 
valid in the n-dimensional case where is an eigenvalue with multiplicity 
fc > 1. In particular, one can still show that, even without condition iii) of 
Theorem II. H the local singular locus of X is still a smooth A;-dimensional 
manifold. However, it is more difficult to show that there is still a local 
regular invariant (n — /c)-dimensional foliation. If one can show the existence 
of this regular regular invariant foliation, then one can drop condition iii) 
from the statement of Theorem I 1 . 1 [ because it is a consequence of the first two 
conditions. Maybe it is possible to use the techniques of Belitskii-Kopanskii 
[1] together with a kind of desingularization of the first integrals in order 
to show the existence of an invariant regular foliation, but we don't have a 
proof so far. 
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